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1.  INTRODUCTION 

Some  work  has  been  done  on  calibrating  the  Large  Scale  Gap  Test 
(LSGT)  for  pentolite  donor  charges.  The  calibration,  which  is 
preliminary,  is  reported  elsewhere’^.  This  report  is  concerned  with 
the  methodology  of  data  acquisition  and  reduction  for  gap  test 
calibration.  One  reason  for  studying  the  methodology  is,  of  course, 
to  obtain  the  best  possible  calibration  of  the  LSGT.  Another 
reason  is  to  find  if  the  methods  are  sufficiently  sensitive  and 
accurate  to  detect  subtle  changes  in  the  velocity  of  propagation  of 
the  shock  in  FMMA  cylinders.  That  is,  does  the  shock  front  progress 
in  a  steady  manner,  or  does  its  velocity  chm'ige  abruptly  at  certain 
distances  along  the  cylinder?  There  are  reasons  for  suspecting  that 
the  acceleration  is  not  a  smooth  function  of  time,  T  (or  of  distamce, 
X),  One  of  the  primary  reasons  for  suspecting  this  is  that  the 
pressure  decreases  from  about  170  kilobars  (kb)  at  the  interface 
between  the  PMMA  and  the  pentolite  to  about  17  kb  at  5  cm  frca  the 
interface.  It  is  improbable  that  this  rapid  attenuation  is  achieved 
by  a  smooth,  continuous  process.  In  the  first  several  cm  of  travel, 
the  shock  front  is  probably  overtaken  by  rarefactions  having 
amplitudes  of  many  kilobars.  These  large  amplitude  waves  should 
produce  abrupt  changes  in  the  shock  velocity  and  its  pressure. 

Kolcky®  has  discussed  wave  reflection  and  reinforcement  behind 
shock  fronts  in  explosive  loaded  cylinders.  With  regard  to  the 
fracture  along  the  axis  of  such  cylinders,  he  notes  that  the  fracture 
does  rK>t  fall  off  uniformly  with  distance  from  the  charge,  but  passes 
through  a  maximum  and  then  decreases  rapidly.  This  is  the  case  for 
a  charge  having  a  smaller  diame^er  than  the  loaded  cylinder.  For 
the  gap  test  the  two  diameters  are  the  same,  so  the  fraccuring 
differs  from  that  reported  by  i'olsky.  In  the  region  within  two  cm 
of  the  charge/PMHA  Interface,  all  of  the  FMMA  is  shattered.  Starting 
at  about  two  cm,  some  of  the  material,  near  the  axis  is  recovered  in 
the  form  of  a  tapered  cylindrical  slug.  For  12.5  cm  long  cylinders 
loaded  with  tetryl,  the  slug  remains  attached  to  the  far  end  of  the 
cylinder^.  The  fractures  meet  the  wall  of  the  cylinder  about  10  cm 
from  the  explosive.  For  pentolite  loading,  the  central  slug  is 
always  broken  loose  from  the  recovered  portion  of  the  cylinder  as  if 
the  tension  waves  had  increased  in  strength.  The  breakout  of  the 
fractures  is  still  about  10  cm  from  the  explosive.  This  breakup  of 
the  cylinders  is  another  reason  for  expectir.g  the  existence  of  relief 
waves  of  fairly  large  amplitude.  Again,  these  relief  waves  should 
produce  observable  effects  on  the  velocity  of  the  shock  front. 

Curvature  of  shock  fronts  in  PKMA  has  also  been  studied*.  These 
indicate  that  the  radius  of  curvature  of  the  front  changes  abruptly 
at  2.5  cm  from  the  interface.  This  is  undoubtedly  produced  by  the 
interaction  of  finite  amplitude  relief  waves  interacting  with  the 
shock  front. 

Most  of  the  previous  calibration  work  at  NuL  was  done  with  a 
35  mm  streak  camera.  We  now  have  a  70  mm  camera  so  that  any  changes 
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in  acceleration  can  be  detected  more  readily.  Hence  the  record.^ 
obtained  in  the  calibration  prograa  for  the  lot  of  pentoiite 
we  have  been  using  have  been  exaained  exhaustively.  This  stuviy  has 
shed  some  light  on  the  subject  of  shock  wave  attenuation  3n  the 
PMMA  cylinders  and  has  indicated  how  the  calibration  could  be 
obtained  more  accurately, 

2.  DESCRIPTIOH  0?  PLOW  Ifi  THE  PMMA  CYLIHDER 

The  attenuation  of  the  shock  ir-ight  be  approxiisated  by  using 
the  theory  of  coopressible  flow.  The  flow  is  two  diEensional,  tiise  t 

dependent,  and  axially  sjn^stric.  Fbrthereore,  the  flow  behind 
the  shock  front  is  subsonic.  Because  the  fr<»it  is  curved,  the 
flow  behind  it  is  rotational  and  non-isentropic.  Bo  anal^ic 
solutions  are  likely  to  be  found  for  such  flow,  even  for  a  satertal 
having  a  sinple  equation  of  state.  The  probles  is  even  core 
explicated  because  the  effects  of  rigidity  rcust  be  included  in  the 
constitutive  relations  for  FKHA. 

The  probi«n  can  be  attacked  by  using  finite  difference  codes, 
and  at  least  two  attcB5)ts  have  been  eade*^?® ,  This  nethod  is  rela¬ 
tively  expensive  even  with  a  modern  digital  coetputer.  Fine  zoning 
is  requii'ed  if  scall  changes  in  the  shock  velocity  are  to  be 
obscirved  in  the  results.  The  attempts  at  solving  the  problem  iised 
relatively  coarse  zoning  so  that  the  results  are  not  definitive. 

One  fau3.t  of  the  codes  is  that  they  cannot  take  into  account 
the  complicated  fracturing  of  the  nCKA  iimecilately  behind  the  siiock 
front,  (Fracturing  is  an  energy  absorbing  process  which  should  be 
taken  into  account  in  the  energy  balance  during  the  exputations , } 

This  fault  makes  it  questionable  if  the  use  of  codes  will  give  the 
desired  results  in  the  foreseeable  future.  Thus  the  problem  remains 
a  challenge  to  both  the  experimenter  and  the  theorist, 

3.  EXEERIMEBrAL 

The  experimental  arrangement  is  the  same  as  that  used  in 
previous  calibrations®.  Briefly,  a  cylindrical  sample  of  FKXA  is 
shocked  by  detonating  a  charge  of  5C/50  pentoiite  (po  *  1.56  g/cc) 
in  contact  with  it.  The  diameter  of  the  charge  Is  5.08  cm  and 
consists  of  two  pellets,  each  2.5^  cm  thick,  and  Is  point  initiated 
so  that  the  detonation  is  axially  symmetric.  The  sample  has  flat 
surfaces  milled  and  polished  on  its  cylindrical  surface  through 
which  the  streak  camera  views  the  light  source.  The  FHHA  is  5,08  cm 
thick  between  the  flat  surfaces^  these  surfaces  are  0.3  to  0,5  cm 
wide  so  that  the  diameter  of  the  cylinder  is  slightly  greater  than 
5.08  cm.  As  the  shock  travels  through  the  sample,  the  light  source 
is  shuttered.  There  results  a  streak  on  the  film  which  relates 
time  and  distance  for  the  shock  front,  A  camera  record  taken  in 
this  way  is  known  as  a  shadowgraph. 

In  what  will  be  referred  to  as  a  ” regular  shot",  the  F5<MA 
sample  is  10  cm  long,  or  longer.  One  fault  of  these  shots  is 
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th*t  they  give  records  uiilch  have  very  little  curvature  so  that 
accurate  di?ferent54ition  is  difricvdt.  Also  for  these  shots^,  the 
objective  len*?  of  the  casera  is  focused  on  a  scale  placed  beside 
the  sa^le  and  at  the  sace  gecBetrlc  distance  frees  the  csseera  as  the 
axis  of  the  charge «  This  results  in  soee  error  because  the  scale 
factor  for  distance  should  be  that  appropriate  for  the  optical 
distazM^e  to  the  center  of  the  I99!A  sa^le>  not  that  for  the  geometric 
distance.  IRiese  coeceients  apply  to  the  regular  ^K>ts  fired  for  this 
work  and  all  previously  reported  work  which  utsd  the  35  *ss  streak 
caeera.  Results  previously  reported  baseo  on  shot  records  froe 
another  70  sas  caeerl?  weiv  obtained  with  proper  focusii^  snd  the 
appropriate  scale  factor?;,  but  not  with  parallel  light, 

Oce  consequence  of  uslzig  P>«<A  ssesple-  10  ca  long  (or  longer) 
is  that  converging  li^t  ssist  be  used  because  the  usable  di«seter 
of  the  objective  lens  is  less  8.8  ce,  the  disneter  of  the 

yiewing  port.  With  convergent  light  the  record  represents  the 
progress  of  different  portions  of  the  3  dissensicsial  shock  front. 
Illvainatlcg  the  event  with  parallel  light  shmxld  give  shadowgraphs 
which  represent  the  progress  of  the  shock  along  the  axle  of  the 
cylinder.  These  conslderatioriS,  along  with  difficulties  in 
differentiating  the  data,  led  to  a  change  in  the  e^t^risental 
arrangecent  so  that  the  c«era  viewed  only  5  cs  of  HQCA.  Thus  one 
shot  covered  the  muss  of  distance,  O  to  5  cs#  aod  another  shot 
covered  5  to  10  cr.  There  is  sose  overlap  between  the  two  shots 
so  that  the  resislts  can  be  Joined  In  a  reasonable  asazuter.  Because 
the  fields  of  view  have  ^Raaller  diaensions  than  the  objective  lens 
of  the  eaaera,  parallel  li^t  can  be  used.  The  caaera  is  focused 
on  a  scale  placed  behind  2.5^  cs  of  mm,  followlz^  which  the 
eaq;>loslve-BtU  assesbly  is  placed  so  that  the  eaaera  is  focused  on 
the  axis  of  syaaetiy.  Thus  the  scale  factor  for  distance  Is  correct 
for  the  optical,  distazace  to  the  center  of  the  sa^Ie  azjd  eaaera 
fOcuaizig  azvi  lighting  are  eptiaua  for  zdiadowgraphy.  These  shots 
are  referred  to  as  "close-up  shots"  in  the  following. 

k.  XUTA  REZXTCTIOK 

^le  eaaera  records  are  digitised  by  using  the  Uz^irersd 
Telereader.  The  output  of  this  s»chlne  is  autoaatlcelly  puncl^i 
into  ZBM  cards  as  tiK  zwoidzer  of  "ecunt&=^  starting  froe  soae  arbl- 
trazUy  selected  origin.  A  coe^iter  progrSB  is  used  to  convert 
counts  into  tiae  and  distance  by  use  of  the  appropriate  reduction 
factors.  These  data  are  then  differentiated  in  order  to  detersine 
the  shock  velocity.  The  fila  record,  data  cards,  snd  sost  of  the 
printed  counter  output  are  stored  for  future  reference. 

The  azrt  of  differentiating  such  data  has  been  discussed  In 
Boae  detsil  In  a  previous  report’.  As  before,  the  data  are  r»t 
evezily  spaced  In  either  variable  so  that  ordlziary  szsoothing  and 
differentiatir^  fozrulas  cazuK>t  be  used.  Slore  cccplicated  formulas 
cazi  be  used,  but  no  computer  code  incorporating  then  was  svailable 
at  the  tise  that  this  werk  was  initiated.  As  in  Beferezjce  7, 
the  data  were  studied  with  the  aid  of  a  spline  function  which  had 
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been  programmed  for  the  IBM  7090.  This  function  is  a  set  of  cubic 
equations,  '•ach  being  fitted  to  a  subset  of  the  data.  Adjacent 
cubic s  are  Joined  so  that  the  overall  function  is  smooth  at  the 
Joints.  The  first  and  second  derivatives  are  also  required  to  be 
smooth  at  the  Joints.  The  computer  program  evaluates  the  parameters 
for  the  spline  function  and  supplies  the  derivatives  at  any  specified 
value  of  the  independent  variable.  When  the  number  of  points  In 
each  subset,  G,  is  3>  the  function  fits  every  point  so  that  minor 
errors  in  the  data  make  the  derivative  irregular.  Larger  values  of 
G  causes  the  program  to  produce  a  smoother  function  so  that  the 
derivative  is  also  smoother.  The  upline  function  appeared 
attractive  because  the  degree  of  smoothing  could  be  controlled  and 
because  it  was  available. 

After  the  data  had  been  differentiated  by  using  the  sp31ne 
function,  it  appeared  that  the  experimental  work  needed  improvement. 
This  led  to  the  "close-up"  shots  mentioned  above.  When  the  cods 
gave  strange  results  for  these  improved  experiments,  the  Gode 
itsslf  became  suspect.  The  spline  function  subroutine  calls  a 
matrix  subroutine  in  order  to  solve  for  the  coefficients.  When 
the  value  of  G  is  small,  3  or  4  in  the  work  referred  to  above,  the 
matrix  is  fairly  large.  Round-off  error  becomes  serious  so  that 
the  results  are  not  reliable.  Replacing  the  single  precision 
matrix  subroutine  with  a  double  precision  subroutine  improved  the 
results,  i,e.,  the  derivative  was  less  erratic.  Because  we  have, 
at  most,  4  significant  figures  in  the  data  from  the  Telereader, 
it  seemed  unwise  to  use  a  numerical  method  which  required  double 
precision.  For  this  reason,  simpler  methods  have  been  used  to 
differentiate  the  data.  These  are  described  in  the  following. 

When  the  data  are  evenly  spaced  in  the  independent  variable, 
simple  formulas  are  available  for  differentiating  numerical  data*. 

In  order  to  reduce  the  effect  of  noise  in  the  results,  other 
formulas  are  available  for  smoothing  the  data  prior  to  computing 
the  derivative  .  These  formulas  have  to  be  applied  with  a  certain 
amount  of  caution,  else  subtle  trends  in  the  X,  T  data  will  be 
suppressed.  In  the  work  reported  here,  the  data  are  not  equally 
spaced  in  either  X  or  T.  Hence  we  must  use  more  complicated  formulas 
for  smoothing  and  differentiating  the  data.  These  formulas,  as 
well  as  those  for  evenly  spaced  data  referred  to  above,  are  based 
on  interpolation  functions  which  are  used  over  a  subset  of  the  data. 
These  functions  are  usually  low  order  polynomials,  the  coefficients 
of  whiGn  are  determined  by  the  use  of  computer  codes.  In  order  to 
preserve  as  much  of  the  character  of  the  data  as  possible,  poly- 
ncsnials  of  the  first,  second  and  third  degree  have  been  used  in 
this  work.  The  number  of  points  in  a  subset  was  either  3  or  5. 
Actually,  formulas  (in  the  usual  sense)  are  not  used.  The  poly¬ 
nomial  is  fitted  to  the  subset  of  data  in  a  least  square  sense  by 
a  computer  subroutine.  This  method  may  not  be  as  efficient  as  one 
based  on  the  use  of  formulas.  It  is,  however,  flexible,  and  requires 
a  minimum  of  computer  programming,  A  self-contained  subprogram  for 
smoothing  and  differentiating  data  was  obtained  from  R.  T.  Nelson,  Jr. 
at  the  Naval  Air  Test  Station,  Patuxent  River.  Maryland  .  Results 
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from  the  use  of  this  method,  es  well  as  from  the  use  of  the  spline 
function  and  a  method  based  on  subdividing  the  data  for  linear 
fits  are  given  in  the  following. 

5.  REGULAR  SHOTS 

A  series  of  four  shots  was  fired  using  identical  procedures  and 
COTiponents.  Plots  of  the  shock  velocity,  U,  as  a  function  of  the 
distance  from  the  explosive-PMMA  interface,  X,  were  obtained  by 
using  the  spline  function,  first  with  G  »  3  and  then  with  G  «  4 
(that  is,  3  and  4  points  per  subset  of  data).  For  Shot  578, 

Fig,  lA,  the  two  cu2nres  nearly  coincide  over  the  whole  range  of  X. 
Both  curves  have  a  "hump”  in  the  neighborhood  of  X  3*0  cm.  The 
two  curves  for  Shot  579  also  agree  reasonably  well,  but  the  hump 
is  not  so  much  in  evidence,  see  Fig.  IB,  Data  for  Shot  580  were 
assumed  to  be  of  poor  quality  because  the  curve  for  G  «  3  has  an 
oscillatory  component,  see  Fig.  1C.  For  this  value  of  G,  the  spline 
function  passes  through  each  point j  see  the  description  of  the  func¬ 
tion  above.  Any  errors  in  the  observations  tend  to  become  more 
prcMulnent  in  the  derivative  of  the  date.  Finally,  the  results  shown 
in  Fig.  ID  imply  that  the  data  for  Shot  581  are  poor  because  the 
curve  for  G  »  3  has  a  large  oscillatory  component.  At  the  time 
these  data  were  processed,  there  was  no  explanation  for  the 
roughness  of  the  U  vs  X  curves  for  the  last  two  shots.  Part  of 
the  trouble  turned  out  to  be  roundoff  error  in  a  matrix  subroutine 
as  mentioned  above.  There  must  be,  however,  some  noise  in  the 
records  for  the  last  two  shots,  particularly  Shot  581,  because 
alternate  treatments  do  not  give  smooth  curves. 

The  effects  of  changing  from  a  single  precision  to  a  double 
precision  matrix  subroutine  are  shown  in  Figs.  2A  and  2B,  which 
should  be  compared  with  Figs.  1C  and  ID,  Noise  is  still  present 
in  the  U  vs  X  curves,  es;  dally  for  Shot  581,  Pig.  2B,  with  G  =  3» 
This  shot  zsust  have  had  some  dei'ect  either  in  the  ccmiponents  or  in 
the  assembled  experiment  Itself.  The  record  was  difficult  to  read, 
largely  because  the  film  was  not  exposed  uniformly.  This  could 
result  fiw.  an  inferior  light  source,  or  because  of  dust  in  the  slit 
of  the  camera.  The  data  used  here  are  from  a  second  and  more 
careful  reading  of  the  film.  The  original  set  of  data  contained 
even  more  noise  than  the  second  set. 

The  difi’erentiation  formulas  described  in  Section  4  have  also 
been  used  to  obtain  the  shock  velocity  from  the  data  for  the  four 
shots.  The  data  for  a  shot  were  first  smoothed  by  fitting  a  subset 
containing  five  points  with  a  second  degree  polynomial,  using  a 
least  squares  criterion.  The  ordinate  for  the  3rd  point  was  then 
calculated  to  give  a  smoothed  value.  The  distance,  X,  was  used 
as  the  independent  variable  so  that  the  coefficient  C  in  the 
relation  T»A  +  BX  +  CX*  woiild  most  likely  be  positive.  If  T 
had  been  used  as  the  Independent  variable,  the  coefficient  of  the 
second  degree  term  would  probably  be  negative.  Note  that  the 
smoothed  value  of  T  is  not  necessarily  at  the  middle  of  an  Interval — 
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it  is  the  third  ordinate  of  the  subset.  Each  point  of  the  data 
is  smoothed  in  turn  by  the  computer  code.  The  choice  of  five 
points  and  a  second  degree  polynomial  is  entirely  arbitrary. 

The  next  step  in  the  process  is  to  go  througli  the  set  of  data 
containing  the  smoothed  ordinates,  Tg,  fitting  three  points  at  a 
time  with  a  quadratic  function.  The  fit  is  a  least  squares  fit, 
and  the  coefficients  are  used  to  compute  the  derivative  at  the 
second  point.  The  choice  of  three  points  and  a  quadratic  function 
is  again  arbitrary.  Results  for  the  four  regular  shots  are  shown 
in  Pigs.  3A  through  3D«  In  these  figures,  the  shock  velocity,  U, 
is  plotted  as  a  function  of  X.  U  is,  of  course,  the  reciprocal 
of  the  derivative,  dT/dX,  as  computed  by  the  scheme  described  above, 
where  X  was  the  independent  variable.  Note  that  in  Fig.  3A,  there 
is  a  rough  section  in  the  curve.  This  is  partly  obliterated  by 
smoothing  the  values  of  U“'  by  using  in  succession  subsets  of  the 
data  set,  X,  11“*'^.  A  second  degree  polynomial  and  five  points  were 
used  in  this  smoothing  process,  resulting  in  a  new  set  of  ordinates, 

Us"^ .  These  smoothed  velocities.  Us,  are  shown  in  Figs.  3A  through 
3D  also.  As  noted  above,  the  rough  section  of  the  curve  for  Shot 
578  at  about  X  »  12  mm  Is  now  a  little  smoother.  Elsewhere,  the 
smoothed  and  unsmoothed  cui^es  are  practically  indistin/^uishable. 

The  same  is  true  for  the  two  curves  in  Fig.  3B  and  3C.  Shot  581, 

Pig*  3D,  has  a  constant  velocity  region  at  X  =  10  ram;  this  is  similar 
to  the  step  in  Fig,  3A.  With  this  method.  Shot  58I  is  about  as 
good  as  any  of  the  others. 

This  local  smoothing  and  differentiation  method  (Isd^  is  the 
preferred  method  at  this  time.  The  codes  are  simple  and  the  results 
are  reasonable.  That  is,  there  is  less  noise  in  the  results,  U  vs  X, 
especially  for  X  <  10  mm,  and  it  gives  better  results  for  X  »  0  than 
those  obtained,  for  example,  from  the  spline  function.  We  do  not 
expect  to  record  the  effects  of  the  reaction  zone  of  the  explosive 
by  the  optical  method  used  for  these  shots.  But  we  should  be  able 
to  obtain  a  velocity  at  the  interface  close  to  that  which  can  be 
calcu3ated  by  using  impedance  matching  and  the  Chapmen- Jouguet 
parameters  for  the  explosives.  This  gives  a  value  of  U  of" about 
0.624  cm/nsec^ .  Ail  of  the  Isd  results.  Pigs.  3A  through  3D  give 
values  of  U  equal  to  or  greater  than  0.56  for  X  =  0,  This  is  grat¬ 
ifying  because  the  results  obtained  from  the  codes  will  not  have 
to  be  changed  gignificantally  for  small  values  of  X  in  order  to  give 
the  correct  Interface  value.  It  is  true  that  some  of  the  curves 
obtained  with  the  spline  function  give  values  of  0.56  cm/nsec,  or  ^ 

greater,  at  or  near  the  interface.  But  too  many  of  then  oscillate 
for  values  of  X  <  10  mm.  The  spline  function  subroutine  could 
probably  be  changed  so  that  the  interface  value  of  U  could  be  used  , 

as  an  input  parameter.  This  would  give  a  boundary  condxtion  on 
the  function,  at,  or  near  the  first  data  point.  However,  roundoff 
error  might  continue  to  be  troublesome  so  that  the  Isd  method  will 
probably  remain  the  preferred  method. 
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6.  CLOSE-UP  SHOTS 

The  problem  of  differentiating  data  from  streak  camera  records 
can  be  better  appreciated  by  examining  the  curves  in  Pig.  4.  These 
show  the  relation  between  the  time  T  and  the  distance  of  shock  travel, 
X,  The  origin  for  curve  A  is  at  the  explosive-PMMA  interface  while 
that  for  curve  B  is  5«08  cm  from  the  interface.  Points  on  the  curves 
represent  the  points  read  on  the  Telereader,  whose  output  is  in 
counts.  For  these  data,  a  count  represents  0.03  nim  in  the  distance 
direction  and  0,007  M-sec  in  the  time  direction.  Points  are  read  at 
spacings  of  about  30  to  4o  counts,  or  0,9  to  1.2  mm,  on  the  film. 

This  gives  40  or  more  points  for  each  of  the  close-up  shots,  see 
Pig.  4.  The  curves  are  nearly  linear,  even  curve  A,  along  which 
the  pressure  drops  from  170  kbar  for  X  «  0  to  about  17  kbar  for 
X  «  5  cm.  This  rapid  attenuation  of  pressure  gives  a  change  of 
slope  which  is  little  more  than  perceptible.  Nevertheless,  these 
data  are  not  as  difficult  to  work  with  as  those  described  above 
where  10  cm  of  shock  travel  was  crowded  into  one  record.  Increasing 
the  size  of  the  plots  referred  to  above  would  help  in  determining 
if  the  derivative  was  continuous  or  not.  Such  a  change  of  scale  can 
be  effected  by  fitting  the  data  with  a  function  linear  in  T  and 
evaluating  the  residuals,  (X  -  Xotserved)*  Sesiduals  are  relatively 
large  for  Shot  715>  Pig*  5A  and  smaller  for  Shot  716,  Pig.  5B,  Each 
small  division  for  the  ordinate  axes  is  0.1  mm  for  Figs.  5A  and  pB  as 
COTipared  to  the  value  of  a  count  which  is  0.03  nui.  The  spatial 
resolution  of  the  streak  camera  is  about  25  lines/mm.  For  these 
shots,  the  magnification  is  about  unity,  so  that  a  count  on  the 
Telereader  is  about  equal  to  the  distance  that  can  be  resolved. 

Thus  the  0,1  mm/di’/lsion  scale  for  the  ordinates  in  Pigs.  5A  and 
5B  is  reasonable  considering  the  resolution  of  the  camera. 

An  interesting  feature  of  Pigs,  5A  and  5B  is  that  sets  of  the 
residuals  can  be  represented  by  straight  lines.  These  straight 
lines  fit  the  points  in  each  set  to  very  nearly  within  ±0.03  nna, 
or  within  one  count.  In  Pig,  5B  the  fits  are  not  as  good,  deviating 
by  something  like  ±  2  counts,  or  ±  0.06  mm.  The  velocity  of  the 
shock  front  can  be  obtained  from  the  coefficients  of  the  original 
linear  relation  between  X  and  T  and  the  linear  relations  between 
the  residuals  and  T,  In  preference  to  such  a  process,  the  plots  are 
used  to  select  sets  of  data  for  linear  "its.  Table  1  shows  the 
results  of  these  fits.  In  the  first  3  columns  are  shown  the  values 
of  T,  X  and  the  residuals.  For  each  set  of  data.  Col,  4  gives  the 
standard  deviation  of  the  linear  fit,  cf.  Col,  5  the  velocity,  U, 
and  Col.  6  the  standard  deviation  of  the  velocity  Uy.  For  the  first 
4  sets  of  data,  the  only  residual  greater  than  0.03  rm  is  in  the 
second  setj  it  is  0.036  mni.  Thus  the  fits  are  as  good  as  one  can 
reasonably  expect,  i.e.  within  ±  1  count  on  the  Telereader.  For 
T  ^  11.73  Jisec  the  fit  is  considerably  worse j  these  data  are  frc«n 
Shot  716,  see  Pig,  5B  for  the  residuals  for  the  overall  linear  fit. 
The  value  of  Of  is  O.O62,  or  2  counts  on  the  Telereader,  and  the 
largest  residual  (absolute  value)  is  0.131>  or  about  4  counts.  This 
implies  that  the  record  of  Shot  716  was  of  lower  quality  than  that 
for  Shot  715,  especially  over  its  first  third. 


7 


NO  UR  70-265 


>o  <  - 

N  5- 


<0 

o 

z 

< 

<o 

</> 

»- 

O 

X 

in 

2 

2 


u. 

O 

i/i 

H* 

Ul 

i/% 

«D 

r) 


i/% 

Ui 

Z 


X 

O 


U  3 
O  X 


<  «• 

o 

X  a 

o 

Otti 

• 

«•«  » 

C' 

O  3- 

>•  u 

K-  w 

mu 
X  m 


i 


_  u 


o 

o 


4«lk 

V)  •» 

.4 

< 

Sx 

N  4  M 

0#4  o  M  B 

090000 

v> 

•  •  4  •  •  • 

iU 

OOOOOO 

flC 

t  1  1 

yj 

u 

c 

B  r»  O  «44  M 

<  X 

BO  C  Ml  •>  M 

n-X 

B  O  O  M  M  O 

«l 

•  ••••• 

0«4  O  B  «»  M 

O 

OOOOOO 

o 

♦ 

o 


o 

o 


0000000000090COPO 


cooooocooccoo 


ooooooooooooooooo 
I  I  I  1  i  t  T 


*4  «0  ^  f# 

•  •  •  •  •  s 


oooooooooocoo 
I  t  t  I  i  I 


^4  t*  t*t  ^t>*r4 


m 


s 

o 


c 

o 


o 

o 


V) 

o 

z 


¥i 

o 


o 

i/> 


D 

U) 

Ul 


•<  — 

5t: 

•4  lb 


«i 

< 

w 

m 


u 

Si 

«c 


t\i 

*«  lA 


o 

o 


o 

e 


p 

o 


»4 

p 

o 


♦ 

#x 

o 


o 

o 


«4^e«Mvo  ooo4»o«*  •^0w04*p  4<%ao<^o  2^222^ 

cie  oooeeo  oooooo  oooooo  opoopo  ooooee  oo««04*oopoqp 

*•  #••••#  ••##•4  #•••••  #•••##  •••### 

eo  eeeoae  eoeooo  eooeoo  oooooo  oooooo  oooooocoooo 

^  ««  T  ■■  ■*  CCaft  Itl  f  lilt 


€4 

•  « 

oo 


o 

r  • 

OP 


I  I 


o 

O  oo  mpo 

«>4««0  OMO 
#4  p«  ^  M  ^  tfl 


|iv  o  o  oo  o 
ooo4-*.r- 

0»  o  <vo  ^> 

!••«€• 

OOP  •'-MO 


.  ^OOO 
f  t  f 


M  Of-  *-••►  f 
O  O  o4  *4  o  O 
<•  Mr«OM  O 
•  ••«»«• 
0  00  4-0^ 
#4  #4  M  M  M  M 


OOOOOO 
I  f  I 


V  •  •  -  -  - 

OOO  O  OO 
I  t  I 


Mo#^<o  Moe*4rM 
OWOO-***  M0#*4«00 

#^MOOM  MOPO^^ 


MMOOOP  OOM  ’•44ir»  npOMO*  f^a«*Or»PM 

oof-oi^«  ^.4#r»oo  «rt44^^«»2#  ^'2 

#«p4Mtflf^O  #>MOOMO 

0  »•  4  9  9  9  9  •  9  9  •  •«•••*  9  9  9  9  9  9  9  9  9  9  9  9 

OOOOOO  «444o4MMM  ^  ^  ^  c^OOOOO  mooooo 


OVWOVVOVO** 

oooooocoooo 
*  II  If 


OOMOOMOC^BO 

BOOOMUMB44MO 


Ortor.POMB^MK 
r»  044mbbbbo.«4B 
4«MMMMMMMinfnn 


; 


)).««4  yr^iio  -o«o«2 

14.021  S«.Si»  -0.004 

14.303  53.205  -0.033 


NOLTR  70-265 


The  results  of  the  process  described  above  are  shown  graphically 
iri  Pig.  6.  So  far,  values  of  U  have  been  calculated  only  over  the 
values  of  X  spanned  by  the  solid  lines.  In  some  cases,  the  shock 
velocity  apparently  changes  discontinuously,  as  between  sets  1  and 
2,  see  Table  1,  and  between  sets  2  and  3»  No  data  were  left  out 
at  these  "joints”.  One  data  point  was  left  out,  however,  between 
sets  3  atid  between  sets  5  and  6,  see  Fig,  5A,  Severa'  points 

were  left  out  between  sets  6  and  7.  There  remains  the  problem  of 
connecting  the  lines  for  constant  values  of  U  in  Fig,  6,  at  least 
where  data  have  been  excluded  from  the  fits.  Over  these  regions, 
the  shock  velocity  may  be  contihuous  while  the  acceleration  under¬ 
goes  large  changes. 

Even  a  casual  reader  must  have  noted  that  the  curves  in  Figs.  5A 
and  5B  are  parabolic  in  shape.  Therefore  they  could  be  fitted  in 
a  least  squares  sense  by 

R  =  A  +  BT  +  CT*, 

Where  R  represents  the  residuals,  (X  -  Xe),  where  Xe  is  an  observed 
value  of  the  position  of  the  shock.  Values  of  X  are  from  the 
linear  fit. 


X  =  a  +  bT 

If  we  assume  the  residuals  are  fitted  exactly,  the  equations  can 
be  combined,  giving 


Xg  a  (a  -  A)  +  Oi-  B)T  c  CT®, 

so  that  Xe  and  T  are  related  exactly  by  a  quadratic.  TMs  method  of 
smoothing  data  is  similar  to  a  method  developed  by  A,  T.  Doodson 
and  discussed  by  Hartree^* .  The  latter  cautions  that  smoothing 
data  may  result  in  errors.  He  states; 

"The  main  purpose  in  carrying  out  a  process  of  smoothliig 
must  therefore  be  to  achieve  smoothness,  not  accuracy. 

The  contexts  in  numerical  analysis  in  which  smoothness 
is  a  prime  requirement  are  not  many,  so  that  such  a 
P2x)ce8s  is  not  often  required.  But  occasionally  it  is 
difficult  to  make  satisfactory  progress  without  one." 

Differentiating  streak  camera  records  is  a  process  in  which  smoothing 
is  required  in  order  to  make  progress. 

Rather  than  evaluate  the  five  coefficients  in  the  equations 
given  above,  the  X,  T  data  were  fitted  directly.  The  best  quadratic 
fit  by  least  squares  is 

X  =  0.36  +  5.33  T  -  0.09T* 

so  that  the  velocity  is 
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u  =  5.33  -  0.18T. 

This  is  an  alternate  interpretation  of  the  record  for  Shot  715. 
A  similar  interpretation  of  Shot  716  gives 

X  =  0.02  +  3.28T  -  0.009T* 


so  that 


u  =  3.28  -  0.018T. 

These  interpretations  of  the  data  for  Shots  715  and  716  require 
the  acceleration  of  the  shock  front  to  change  by  a  factor  of  10 
after  traveling  50  imn.  One  has  the  choice  of  believing  that  there 
is  this  one  big  change  in  the  acceleration  (at  X  =  50  mm),  or  that 
there  are  several  smaller  changes  (as  in  Pig.  6)  or  thaftne  ?ccel- 
eration  is  smooth  as  has  heretofore  been  assumed.  This  is  not  an 
enviable  situatioa. 

The  situation  described  abov^  can  be  complicated  even  further 
by  studying  the  residuals  from  a  quadratic  fit.  This  leads  one  to 
believe  that,  for  example,  a  cubic  should  be  used  to  fit  the  data 
for  Shot  715.  Thus  one  is  led  alor.g  the  path  to  polynomials  of 
increasing  degree.  Because  there  are  no  physical  reasons  for  ejqpecting 
the  shock  path  to  be  described  by  higher  degree  polynomials,  that 
route  has  not  been  followed. 

7.  LOCAL  SMOOTHING  AND  DIFFERENTIATION  OP  ‘’CLOSE-UP"  SHOTS 

Data  from  these  shots  have  been  studied  by  using  the  Naval  Air 
Test  Station  subroutine  "Grout”  (see  Appendix  A).  The  data  have 
been  smoothed  and  differentiated  using  first,  second  and  third 
degree  polynomials.  Only  part  of  the  results  are  reported  heres 
those  given  in  Pigs,  JA,  7B,  8A  and  8B  are  based  on  second  degree 

smoothing  over  5  points.  The  smoothed  data,  X,  Tc,  were  then 

differentiated  by  using  the  second  degree  polynomial  over  3  points; 
thus  the  three  points  are  fitted  exactly  and  the  de:.'ivative  is 
evaluated  for  the  second  value  of  X  of  the  3  values  being  used. 

Figures  7A  and  8A  show  results  obtained  in  this  way  for  Shots  715 

and  716  respectively.  Figures  7B  and  8b  show  results  of  smoothing 
the  values  of  U  by  using  a  second  degree  polynomial  over  5  points. 

The  discontinuous  X,  U  values  of  Pig.  6  are  superimposed  on 
Figs.  7A  through  8B.  There  is  some  correlation  between  the  results 
from  the  Isd  process  and  the  discontinuous  results,  see  Figs,  7A 
and  7B.  Smoothing  the  X,  U  data  makes  the  correlation  less 
noticeable.  Here  again  is  the  dilemma  which  is  always  encountered 
when  numerical  differentiation  is  required;  how  much  should  the 
data  be  smoothed?  If  one  is  biased  in  favor  of  a  smooth  shock  path 
(X,  T  curve)  and  smooth  acceleration,  he  will  prefer  the  smooth 
values  of  U  as  given  in  Figs,  'fB  and  SB.  Bias  toward  the  more 
chaotic  condition  as  suggested  by  Figs.  5A  and  5B  makes  one  prefer 
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the  unsffioothed  results.  Pigs,  7A  and  8A,  A  wise  choice  obviously 
depends  on  how  Violl  the  results  can  be  repeated,  which,  of  course, 
depends  on  the  care  with  which  the  experiments  are  performed,  and 
on  the  precision  of  the  camera  and  film  reader.  These  points  will 
be  discussed  in  later  sections. 

For  small  values  of  X,  the  values  of  U  are  off  scale  in  Pigs, 
7A  and  7B«  The  first  several  values  of  X,  U  and  Ug  are  given  in 
Table  2,  The  smoothing  operation  did  not  change  the  value  of  U 


Table  2 


Values  of  Shock  Velocities,  U  and  Ug 
for  Small  Values  of  X  (Shot  715) 


X(inm) 

U(mm/usec) 

Ug  {wm/w 

0.00 

7.12 

7.23 

0.046 

6.64 

6,49 

1.11 

5.88 

5.88 

1.69 

5.50 

5.62 

2.93 

5.51 

5.51 

4.25 

5.55 

5.54 

very  much  over  this  range.  For  the  first  two  points,  both  U  and 
Ug  are  greater  than  the  Chapman-Jouguet  value,  6,24  mm/>tsec  (see 
Section  5) •  The  more  or  less  flat  portion  of  the  curve  where  U  ^ 
5.5  nm/iisec  coiild  possibly  be  the  effect  of  the  reaction  zone. 

There  ere  good  reasons  to  be  doubtful  of  these  results  for  small 
values  of  X:  (a)  the  numerical  methods  cannot  give  accurate  results 
at  the  end  of  the  set  of  data  auid  (b)  the"flats"  on  the  PKMA  samples 
Y^re  not  good  neaj*  their  ends.  These  samples  should  be  of  optical 
quality;  this  is  difficiilt  to  achieve  in  a  non-optical  shop.  The 
flats  were  polished  by  hand  which  introduced  curvature,  especially 
near  the  ends  of  the  flats. 

For  the  same  reasons  as  given  above,  the  values  of  U  for 
X  50  am  In  Figs.  7A,  7B,  8A,  and  8b  are  probably  not  reliable, 

8.  REPROD»;CIBILITY  CF  U(X)  FROK  THE  POUR  REGULAR  SHOTS 

The  methods  which  have  been  used  for  obtaining  U(X)  from  the 
data  are  objective.  That  is,  they  are  numerical  schemes  which  give 
the  same  results  irtien  a  given  set  of  data  are  treated  a  second 
time.  Admittedly  the  "segmented  fits”  to  the  data  of  Shots  715 
and  716  were  based,  to  some  degree,  on  a  subjective  selection  of  the 
subsets  of  data.  The  other  methods  were  objective,  and  of  these 
the  local  smoothing  and  differentiation  method  has  been  chosen  for 
a  study  of  reproducibility.  This  choice  is  also  admitted  to  be 
subjective.  If  we  had  enough  Information  to  permit  an  objective 
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selection  of  a  method  for  reducing  the  data,  vre  would  undoubtedly 
have  all  the  information  needed  for  the  calibration  of  the  LSQT 
with  no  i\irther  work.  In  the  following  discussion,  reproducibility 
is  based  on  the  results  from  the  Isd  scheme,  using  quadratic  fits 
to  5  points  for  smoothing,  and  to  3  points  for  differentiation. 

This  choice  is  arbitrary:  it  is  rationalized  on  the  basis  that  it  is 
simple,  it  smooths  the  data  a  minimum  amount,  and  the  results  appear 
to  be  reasonable  (an  admitted  subjective  appraisal). 

Results  for  the  four  shots  discussed  previously  are  shown  in 
Fig.  9  as  cuir/es  relating  U  and  X.  A  value  of  U  is  plotted  for  each 
pair  of  values  of  X  and  T  in  the  original  sets  of  data.  For  each 
shot,  the  points  are  connected  by  straight  linesj  this  is  the  way 
the  CalccMnp  Plotter  produces  "curves”.  While  comparing  these  curves 
it  must  be  remembered  that  the  data  reduction  technique  may  give 
unrepresentative  results  at  each  end  of  a  set  of  data.  This  is 
inherent  in  numerical  methods  for  interpolatioris,  smoothing  and 
differentiation.  Hence  the  divergence  of  the  results  for  the  four 
shots  at  X  *  0  is  not  significant.  Otherwise,  divergence  of  the 
curves  represents  lack  of  reproducibility  from  shot  to  shot.  At 
X  =  5  mm,  the  total  spread  n  U  is  0.l4  mm/»ir>ec,  the  average  value 
of  U  is  about  5.35,  so  the  spread  is  ±  1.35^.  At  X  =  25  mm,  TJ  =  4.29, 
and  the  spread  in  values  of  U  are  ±  0.11  mm/psec,  or  ±  2.6^. 

When  tetryl  was  lest  calibrated  in  the  LSGT*  ,  error  in  U  at 
X  «  5  mm  was  estimated  to  be  about  55^  in  the  average  of  five  shots. 
For  U  “  5.3  mm/psec,  (X  approximately  5  mm)  the  range  in  values  of 
U  was  estimated  to  be  about  ±  0.13  or  ±  2.5S^.  Hence  the  spread  in 
the  shock  velocity  for  the  four  shots  discussed  above  is  not 
unusually  large.  In  fact,  the  reproducibility  of  the  results  frwn 
these  four  shots  is  entirely  ccxaparable  to  that  obtained  in  previous 
calibration  work.  Of  course,  such  reproducibility  is  the  c (Elective 
effect  of  (a)  replication  of  booster  pellets,  (b)  adequacy  of  optics 
and  camerc^  and  (c)  methods  of  data  reduction. 

Even  though  the  reproducibility  is  good  for  the  four  regular 
shots,  the  magnitude  of  the  spread  of  the  values  of  U  make  it  impos¬ 
sible  to  detect  any  subtleties  which  could  be  the  result  of  wave 
interaction.  That  is,  the  breaks  in  the  U  vs  X  relation  shown  in 
Fig,  6  are  not  detectable  within  the  spread  of  the  results  of  the 
four  regular  shots.  Any  argument  about  wave  interaction  must  be 
based  on  the  results  from  the  close-up  experiments,  .Shots  715  and 

716. 


9.1  IMPROVING  THE  SHOCK  VELOCITY  METHOD  OP  CALIER.ATIKG  GAP  TESTS 

Some  comments  have  been  made  in  the  preceding  text  about 
improving  the  calibration  of  the  gap  test.  These  improvements  are 
discussed  in  r-jre  detail  in  the  following  two  sections.  When  it 
is  feasible,  these  improvements  will  be  included  in  future  work. 
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9.2  IMPROVEMEWrS  IN  COMPONENTS 

A  significant  improvement  can  probably  be  made  by  more  careful 
preparation  of  the  flat  surfaces  on  the  PMMA  cylinders.  In  the 
past,  these  have  been  prepared  by  machining  and  hand  polishing  in 
an  ordinary  machine  shop.  Such  a  shop  cannot  be  expected  to  produce 
specimens  of  good  optical  quality.  The  hand  polishing  is  probably 
the  most  difficult  operation.  There  is  a  tendency  to  produce  a 
curved  surface,  the  radins  of  curvature  of  which  is,  of  course, 
greater  than  that  of  the  cylinder.  Thus  the  two  flats  become  the 
equivalent  of  a  barrel  lens.  Near  the  ends  of  the  sample,  the 
surfaces  become  curved  toward  the  axis  of  the  charge  also.  That  is, 
the  two  surfaces  are  closer  together  at  the  ends  of  the  cylinder. 
Optical  distortion  because  of  this  defect  may  have  been  the  reason 
for  stating  that  measurements  could  not  be  made  for  X  <  0.5  cm  . 

For  the  next  calibration,  PW!A  specimens  will  be  used  which 
were  prepared  in  the  Optical  Shop  of  the  Bureau  of  Standards.  With 
these,  it  is  hoped  that  data  from  near  the  explosive/PMMA  interface 
will  be  mere  reliable. 

In  the  past,  exploding  bridge  wires  were  used  as  the  sources 
of  light.  Inexpensive  lenses  were  used  to  provide  nearly  parallel 
light  through  the  sample  because  they  are  destroyed  by  the  blast. 
Better  lenses  should  be  used;  someidiat  better  lenses  have  been 
obtained  for  work  in  the  immediate  future. 

Even  better  lighting  should  be  used  for  these  experiments.  Some 
modification  in  the  bomb-proof  would  permit  placing  the  light  source 
outside  the  bomb-proof.  This  would  permit  the  use  of  a  good  lens 
so  that  practically  parallel  light  would  be  obtainable.  A  laser 
could  eventually  be  incorporated  into  the  system  for  even  better 
light. 


The  optical  system  can  be  improved  by  usir.g  better  material  in 
the  port  in  fi-ont  of  the  camera,  FMMA  is  usually  used  in  camera 
ports  for  this  kind  of  work  at  NOL.  Many  installations  use  gloss, 
and  exploratory  work  has  shown  us  that  glass  in  the  port  upgrades 
the  quality  of  the  image  in  the  streak  camera.  Glass  will  be  used 
for  future  calibration  work. 

One  "component"  over  irtiich  we  have  no  control  is  the  Telereader 
which  is  used  to  digitiiie  the  data.  This  machine  is  In  another 
Department.  Its  performance  has  been  upgraded  since  it  was  last 
discussed’’'.  This  was  done  by  replacing  vacuum  tube  circuits  by 
solid  state  circuits,  resulting  in  more  stable  operation.  The  optics 
of  the  instrument  are  still  deficient  in  that  the  lighting  is  uneven. 
This  can  cause  trouble  because  the  operator  may  compensate,  for 
example,  by  moving  the  cross  wires  further  into  the  streak  in  brightly 
lighted  axaas  of  the  field  of  view.  Reading  a  record  a  second  time 
on  this  machine  may  give  sigiiificantally  different  results  when 
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the  digital  data  are  processed.  Placement  of  the  wires  on  a  streak 
is  a  subjective  operation,  of  course.  It  should  be  made  less 
difficult  by  improving  the  lighting  in  this  machine. 

k  second  factor  outside  our  control  is  the  quality  and  uniformity 
of  the  pentolite  pellets  supplied  by  the  Naval  Anmninition  Depot, 

Crane,  Indiana.  However,  when  we  feel  that  we  have  reached  the 
practical  limits  of  improvements  in  the  other  components  and  techniques, 
a  special  experimental  lot  of  pellets  can  be  prepared  for  additional 
studies . 

It  has  been  impliclty  assumed  in  the  above  discussion  that  the 
PUMA  is  reproducible  as  a  material. 

9.?  IMPROVEMSKI  IN  TECHNIQUES 

Some  improvements  in  techniques  have  been  mentioned  in  the 
discussion  about  the  close-up  shots.  When  the  field  of  view  is 
restricted  to  less  than  8.9  cm  (3.5  inches),  the  light  through  the 
sample  can  be  parallel.  That  is,  the  field  of  view  is  smaller  than 
the  port,  which  is  8.9  cm  in  diameter.  The  objective  lens  has  a 
focal  length  of  3O.5  cm  (12.0  Inches)  and  a  diameter  of  about  10.9  cm 
(4.3  Inches),  so  it  is  not  limiting  the  field  of  view.  Parallel 
light  in  shadowgraphy  produces  more  “legible”  pictures^ It  also 
instates  that  the  backlight  is  shuttered  by  that  part  of  the  shock 
wave  in  a  plane  which  contains  the  axis  of  the  cylinder,  and  idiich 
is  perpendicular  to  the  optic  axis.  For  converging  light,  as  was 
used  for  the  regular  shots,  the  “shuttering  plane”  at  first  lies 
behind  the  axis  of  the  specimen,  moves  forward  tuitil  it  is  In  frcmt 
of  the  axis  (assiming  that  the  cylinder  is  located  symnetrically 
with  respect  to  the  optic  axis).  Movement  of  the  shuttering  plane 
produces  an  error  in  the  camera  record  wliich  cannot  be  elisdnated 
during  data  processing. 

In  order  to  obtain  superior  records,  the  streak  camera’s 
objective  must  be  focused  on  the  plane  aensioned  above;  the  plane 
containing  the  axis  of  the  cylinder.  This  can  be  done  by  first 
focixslng  on  the  front  flat,  and  then  or;  the  rear  flat,  and  taking 
the  average  position.  Or  a  good  piece  of  FKMA  half  as  thick  as  the 
distance  between  the  flats  cn  the  cylinder  can  be  put  into  the  opticail 
path  while  Uie  camera  is  being  focused. 

Many  details  uucf  be  kept  in  mind  when  an  experiment  is  assembled. 
Obvlottsly  the  optics  must  be  clean.  The  slit  must  be  properly 
aligned  and  adjusted  as  well  as  cleaned.  Light  source,  subject, 
and  camera  should  be  aligned  to  optical  bench  precision  —  even 
though  temporary  stands  are  used  in  the  bomb-proof.  There  are 
sufficient  minor  problems  in  these  experiments  to  cause  ona  to 
despair  of  ever  obtaining  the  ultimate  accxn>acy.  Adding  these  to 
major  problems,  such  as  discussed  above,  makes  it  obvioiis  that  cali¬ 
bration  work  is  a  continuing  project.  Aside  from  the  necessity  of 
Increasing  reflneoont  of  the  experli^nt,  each  new  lot  cf  booster 
pellets  must  be  tested. 
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10.  CONCLUSIONS 

Several  iaprovesents  In  the  data  acquisition  phase  of  gap 
calibration  have  been  suggested.  These  include  the  use  of  better 
optical  conponents  and  better  procedures- 

The  data  reduction  problem  has  been  examined  extensively.  It 
was  found  that  the  spline  function  is  not  useful  in  differentiatiiig 
data  from  the  gap  test  calibration.  This  is  contrary  to  a  previously 
published  opinion  cf  the  author.  Other  Implegentations  of  the  spline 
function  could  perhaps  be  useful. 

Smoothing  and  differentiating  by  numerical  operations  on  subsets 
of  the  data  appear  to  be  the  best  method  of  obtaining  the  slx>ch 
velocity.  These  numerical  operations  eure  somewhat  co^licated  by  the 
fact  that  the  data  are  not  ^ually  spaced.  The  use  of  a  modem 
digital  cceputer  makes  \hit^  complication  trivial  once  the  codes  have 
been  perfected.  The  flexibility  of  this  method  is  a  most  desij^bie 
attribute;  the  aiscunt  of  smoothing  can  be  c<w)trolled  so  -hat  fairly 
subtle  trends  in  the  data  can  be  detected.  On  the  other  rmnii,  the 
data  can  be  smoothed  to  any  degree  desired,  making  it  possible  to 
suppress  any  suspected  noise. 

It  has  not  been  possible  to  give  a  clear  cut  demonstration  that 
the  shock  velocity  vs  X  curve  has  disc<»itinuities  or  that  its  slope 
is  discontinuous.  The  data  for  the  four  regular  shots  do  suggest 
that  the  “bulge"  reported  previously  does  exist.  CoBqpatations  tend 
to  confirm  this  <*servafcion.  Restate  from  Shot  715,  in  which  only 
5.0  cm  of  the  BQiA  cylinder  was  observed,  suggests  that  the  V  » 

U(X)  relation  Is  even  more  complicated.  This  observation  is  being 
checked  during  the  calibration  of  the  next  lot  of  pentolite  pellets. 
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FIG.  lA  SHOCK  VELOCi 
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FIG.  1C  SHOCK  VELOCITY  AS  DETERMINED  BY  A  SPLINE  FUNCTION  FOR  SHOT  580 
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METHOD  SHOT  578 
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METHOD,  SHOTSeO 
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TIME  VS  DISTANCE  FOR  SHOCK  IN  A  CYLINDER  OF  PMMA,  SHOTS  715  AND  716 
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RESIDUALS  ROM  LINEAR  FIT  TO  X,T  DATA  SHOT  716 
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DISTANCE  (mm) 

VELOCITY  AS  COMPUTED  BY  FIHING  SUBSETS  OF  DATA  WITH  STRAIGHT  LINES,  SHOTS  715  AND  716 


SHOCK  VIIOCITY  (inm/»SEC) 
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DISTANCE  (mm) 


FIG.  7A  SHOCK  VELOCITY  AS  COMWTEO  BY  LOCAL  SMOOTHING  AND  DIFFERENTIATING,  SHOT  715 
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DISTANCE  inn) 


m.  n  SHOCK  VEIOCITY  OF  FIGURE  7A  SMOOTHED  SY  USING  A  QUADRATIC  OVER  FfVE  POINTS 


VCLOCI^  <«m/M  SK)  SHOCK  VELOCITY  (m<ii//»sec) 
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FIG.  8A  ^  COMPUTED  lY  LOCAL  SMOOTHING  AND 

differentiating,  shot  7W‘ 


3t>  40 

DISTANCE  (mm) 


FIG.  M  SHOCK  VELOGTY  OF  FIGURE  8A  SA400THED  BY  USING  A  QUADRATIC 
OV«  FIVE  POINTS,  SHOT  714  QUADRATIC 
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FIG.  9  SHOCK 
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Appendix  A 


SUBROUTINE  ?0R  SMOOTHING  AND  DIFFERENTIATING 
UNEVENLY  SPACED  DATA 

A  subroutine  has  been  acquired  from  Patuxent  River  Naval  Air 
Station  which  smooths  and  differentiates  unevenly  spaced  data.  The 
program  is  known  by  the  name  ’’Crout”  and  is  described  in  the  attached 
subroutine  abstract  by  Richard  T,  Nelson>  Jr.  The  description  is 
someirtiat  misleading  with  regard  to  obtaining  the  derivative  of  the 
smoothed  input  data.  If  the  values  of  the  derivative  of  the  smoothed 
data  are  wanted,  the  smoothed  ordinates  must  be  used  in  a  second  call 
of  the  subroutine.  Some  changes  were  made  in  the  code  so  that  it 
w«Kild  run  on  the  IBM  7090  as  a  Fortran  IV  code.  A  listing  of  the 
revised  code  can  be  acquired  by  a  request  addressed  to  the  author 
of  this  report. 

Changes  in  the  Subroutine 

1.  The  variable  FR  has  been  removed  from  the  call.  It  is 
set  to  1.0  in  the  subroutine.  It  could  be  elminated. 

2.  The  common  statement  has  been  deleted. 

3«  The  array  S  is  now  in  the  call,  and  is  dimensioned  £(M). 

4.  All  print  statements  have  been  repr.aced  by  write  statements. 

5.  Tests  on  £(I)  between  statements  15  and  l6,  and  between  30 
and  30  prevent  attempting  to  compute  0**0  in  the  stat<WB^ts  for 
computing  DK(I). 

6.  Following  statement  20,  the  value  of  I  is  saved  for  use 
at  statmunt  50. 
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5  August  1968 
CROUT 


PURPOSE: 

To  sffiooth  M  data  points,  which  may  be  unequally  spaced,  by  the 
method  of  least  squares  and  moving  arc  polynomial »  Differences 
between  the  original  and  smoothed  points  and  the  first  derivative 
of  the  smoothed  points  are  computed. 

RESTRICTIONS: 

See  calling  sequence. 

CALLING  SEQUENCE: 

Call  Grout  (M,  ND,  NS,  FR,  P,  3,  X,  DX,  DIP,  lERR) 

M  Is  the  total  number  of  data  points 
“*  M  %  3>  M  at  NS 

ND  is  the  degree  of  fit 
ND  ^  19,  NIK  NS 

NS  Is  the  number  of  points  to  idilch  the  fit  Is  to  be  made. 

3  <  N3  <  21  and  must  be  odd. 

FR  Is  a  conversion  factor  such  that  the  first  derivative  will 
Be  computed  In  the  desired  terms.  For  example,  for  radar, 
idiere  P  is  In  seconds  and  S  Is  in  feet,  FR  must  be  1  to 
obtain  DK  In  ft/sec.  For  theodolite,  where  F  Is  frame  number 
and  S  Is  In  feet,  FR  must  be  the  number  of  frames  per  second 
to  obtain  DK  in  ft/sec. 

FR  >  0 

P  Is  the  N  dimensional  array  of  abscissas. 

£  Is  the  N  dimensional  array  of  ordinates. 

X  Is  the  M  dimensional  array  of  smoothed  ordinates. 

DX  Is  the  K  dimensional  array  of  the  derivative  of  the  smoothed 
ordinates. 

DIP  Is  the  M  dimensional  array  of  the  differences  between  the 
smoothed,  ordinates  and  the  original  ordinates. 

lERR  will  be  zero  If  no  errors  are  found  In  M,  NS,  ND  and  FR. 

ERROR  RETURNS: 

Grout  -  Totia  No.  of  points  (M)  Is  less  than  3. 
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Grout  -  Total  No.  of  points  (M)  is  less  than  the  No,  of  points 
to  be  fitted  (NS), 

Grout  -  No* of  points  to  be  fitted  (NS)  is  even. 

Grout  -  No*  of  points  f)  be  fitted  (NS)  is  greater  than  21 
or  less  than  3. 

Grout  -  Degree  of  fit  (ND)  is  greater  than  the  No,  of  points 
to  be  fitted  (NS). 

Grout  -  Degree  of  fit  (KD)  is  greater  than  19* 

Grout  -  Conversion  factor  (PR)  is  tero, 

METHOD ; 

When  a  point  is  under  consideration  in  the  sBioothi!ag  operation, 
a  polynomial  of  degree  ND  is  found  for  the  point  and  2K  other 
points  (where  NS  »  2N+1),  Except  for  the  end  points,  (the  first  K 
and  the  last  N),  the  2N  points  consist  of  the  N  immediately 
preceding  and  the  N  immediately  following  the  point  under  consider¬ 
ation,  To  retain  some  significance,  the  abscissa  of  the  point  under 
consideration  is  subtracted  from  Itself  and  the  2N  other  points 
giving  a  new  set  of  abscissas.  An  ND  by  KD  +  1  Matrix  (representing 
the  set  of  NT  -f  1  normal  equation)  is  found  using  the  principle 
of  least  squares.  Prom  the  Matrix  of  normal  equations,  anot!y?r 
Matrix  is  derived  and  solved  by  Grout's  method.  Solution  of  the 
derived  Matrix  yields  the  coefficients  of  the  ND  th  degree  poly¬ 
nomial  which  best  fits  the  point  in  consideration  and  the  2H  other 
points  as  described  above.  Substitution  of  the  new  abscissa  for 
the  point  in  consideration  (which  is  always  zero)  will  yield  the 
smoothed  ordinate  value.  The  smoothed  ordinate  value  for  the  end 
points  (1st  N  points  and  last  N  points)  are  found  by  substituting 
the  new  abscissa  for  the  end  points  Into  the  polynomial  idiich 
describes  the  1st  NS  points  and  the  last  NS  points  respectively. 

The  first  derivative  is  found  by  differentiating  the  polynomial. 

REFERENCE: 

Nielsen,  K.  L.,  "Methods  in  Numerical  Analysis”,  The  Macmillan 

Company  1956 

LANGUAGE: 

This  subroutine  is  coded  in  FORTRAN 
AUTHOR: 

Richard  T.  Nelson,  Jr. 
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